Collective modes in a system with two spin-density waves: the 'Ribault' phase of 

quasi-one-dimensional organic conductors 
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We study the long-wavelength collective modes in the magnetic-field-induced spin-density-wave 
(FISDW) phases experimentally observed in organic conductors of the Bechgaard salts family, fo- 
cusing on phases that exhibit a sign reversal of the quantum Hall effect (Ribault anomaly) . We have 
recently proposed that two SDW's coexist in the Ribault phase, as a result of umklapp processes. 
When the latter are strong enough, the two SDW's become circularly polarized (helicoidal SDW's). 
In this paper, we study the collective modes which result from the presence of two SDW's. We 
find two Goldstone modes, an out-of-phase sliding mode and an in-phase spin-wave mode, and two 
gapped modes. The sliding Goldstone mode carries only a fraction of the total optical spectral 
weight, which is determined by the ratio of the amplitude of the two SDW's. In the helicoidal 
phase, all the spectral weight is pushed up above the SDW gap. We also point out similarities with 
phase modes in two-band, bilayer, or d + id' superconductors. We expect our conclusions to hold 
for generic two-SDW systems. 

PACS Numbers: 72.15NJ, 73.40Hm, 75.30Fv 



I. INTRODUCTION 



In electron systems with broken symmetries, such as 
superconductors or DW systems, quasi-particle excita- 
tions are often gapped, and the only low-lying excitations 
are collective modes. The latter thus play a crucial role 
in various low-energy properties. 

In an incommensurate SDW system, there are two 
(gapless) Goldstone modes: a sliding mode and a spin- 
wave mode, which result from the spontaneous breaking 
of translation symmetry in real ^Bflce and rotation sym- 
metry in spin space, respectively.BH Contrary to the case 
of superconductors, collective modes in DW systems di- 
rectly couple to external probes and therefore show up 
in various experiments. For instance, the sliding mode, 
which is pinned by impurities in real systems, can be de- 
pinned by a strong electric field. This -Leads to non-linear 
conduction, observed in DW systems.El 

The aim of this article is to study long-wavelength 
collective modes in a quasi-lD system where the low- 
temperature phase exhibits two SDW's. The presence 
of two SDW's gives rise to a rich structure of col- 
lective modes, which in principle can be observed in 
experiments.l3 

Our results are based on a particular case: 
the magnetic-field-induced spin-density-wave (FISDW) 
phaseSppt the organic conductors of the Bechgaard salt 
family.QLl These FISDW phases share common features 
with standard SDW phases, but also exhibit remarkable 
properties like the quantization of the Hall effect. We 
have shown that umklapp processes may lead in these 
systems to the coexistence of two SDW's with compara- 
ble amplitudes, which provides a possible explanation of 
the sign reversal of the quantum Hall effect (QHE) (the 



so-called Ribault anomalyHu) observed in these conduc- 
tors. There are two motivations for studying this particu- 
lar case, i) The Bechgaard salts, as a possible candidate 
for a two-SDW system, present their own interest, ii) 
A conductor with two SDW's is in ^neral not easy to 
analyze, even at the mean- field level. Ej The analysis sim- 
plifies when a strong magnetic field quantizes the electron 
motion. 

Nevertheless, we expect our conclusions to be quite 
general and to apply (at least qualitatively) to other sys- 
tems with two SDW's. This hope is strongly supported 
by the similaritiesci that exist between collective modes. 
in two-|SpW conductors and phase modes in two-band ,EJ 
bilayerli3 or d -I- id' (Ref. O) superconductors, and, to a 
lesser extent, pjlasmon modes in semiconductor double- 
well structures.Ej These similarities suggest that collec- 
tive modes in two-component systems present generic fea- 
tures that do not depend on the particular case consid- 
ered. 



A. Umklapp processes in quasi-lD SDW systems 



Consider a quasi-lD conductor with a SDW ground 
state. In presence of umklapp processes transferring mo- 
mentum K (K being a vector of the reciproqual space) , 
spin fluctuations at wave vectors Q and K — Q are cou- 
pled. Thus, the formation of a SDW at wave vector Qi 
will automatically be accompanied by the formation of a 
second SDW at wave vector Q2 = K — Qi , provided that 
Qi 7^ Q2- The case Qi = Q2 = K/2 corresponds to 
a (single) commensurate SDW. Umklapp processes pin 
the SDW whose position with respect to the underlying 
crystal lattice becomes fixed: the sliding mode is gapped. 



For two incommensurate SDW's (Qi ^ Q2)Ji3 the total 
spin-density modulation can then be written as 



(S(r))= Y, S.cos(Q, -r- 



(1.1) 



where r = (na, m\)) (with n, m integers) denotes the po- 
sition in real space (a and h being the lattice spacings 
along and across the conducting chains). 

Even for Qi ^ Q2, the distinction between the two 
SDW's may appear somewhat unjustified since one can- 
not distinguish between cos(Qi • r) and cos(Q2 • r) when 
r is taken as a discrete variable. However, umklapp pro- 
cesses do lead to the presence of two non- vanishing order 
parameters, (cjj|Ck-i-Qij,) and (cj^fCk+Qai), i^ the SDW 
phaseJij This doubles the number of degrees of freedom 
of the SDW condensate, which yields for instance twice 
as many collective modes (as compared to the case with 
a single SDW). Thus, it is natural to speak of two SDW's 
in the ground state of the system. Furthermore, we note 
that cos(Qi -r-l-ffi) and cos(Q2 •r-|-02) are indistinguish- 
able only if d\ — ~Qi (with, again, r being a discrete 
variable). This is precisely the equilibrium condition ob- 
tained by minimizing the mean-field condensation energy 
(see Sec. III). Condensate fluctuations do not in general 
satisfy the condition Q\ = —02, and it is then more ap- 
propriate to view these fluctuations as originating from 
two different SDW's. 

In a quasi- ID system with a single SDW, the wave vec- 
tor Q of the spin-density modulation is determined by 
the nesting properties of the Fermi surface: E(k -|- Q) ~ 
— i?(k), where i?(k) is the energy with respect to the 
Fermi level. [For a perfectly nested Fermi surface, one 
would have E{\i -I- Q) = —E(k).] Umklapp processes are 
important only if both Qi and Q2 = K — Qi are good 
nesting vectors. Otherwise, one of the two SDW's has a 
very small amplitude and can be ignored for any practical 
purpose. 

For a 2D (or 3D) conductor, the geometry of the Fermi 
surface appears to be crucial. Consider the following 
dispersion law, which is linearized in the vicinity of the 
Fermi level: 

Ea{kx, ky) — vp{akx — kp) — 2thC0s{kyb + an) + ■ ■ ■ 

(1.2) 

where k^ and ky are the electron momenta along and 
across the conducting chains, a = -I- (— ) corresponds to 
right (left) movers with momenta close to akp. vp and 
kp are the Fermi velocity and momentum for the motion 
along the chains, ti, the interchain transfe r int egral, and 
b the interchain spacing. The ellipses in (1.2) represent 
small corrections that generate deviations from perfect 
nesting, k is a parameter which parameterizes the shape 
of the Fermi surface. 

Most calculations on quasi- ID SDW systems assume 
K = 0, which corresponds to the Fermi surface shown 
in Fig. ||a. There are two 'best' nesting vectors: Qi ~ 



{2kp,Tr/b — S) and Q2 — {2kp, —ir/b + S), where the small 
correctipHA(|(5| ^ n/b) is due to deviations from perfect 
nesting.tZHla Since Q2 = (4A:_f, 0) — Qi, these two vectors 
are coupled by umklapp scattering if the system is half- 
filled {4kp — 27r/a). Two SDW's with equal amplitudes 
will form simultaneously at low temperature. 

Consider now the case k = n/A in a half-filled band, 
which corresponds to the asymmetric Fermi surface 
shown in Fig. |l|b. This Fermi surface has been pro- 
posed as a good approximaiinn. to the actual Fermi sur- 
face of the Bechgaard salts.Eauj The best nesting vector 
Qi = {2kp,Tr/2b — S) is now non-degenerate. By umk- 
lapp scattering, Qi couples to Q2 = {2kp,—TT/2b + S), 
which is not a good nesting vector. At low temperature, 
two SDW's will form simultaneously, but the one with 
wave vector Q2 will have a vanishingly small amplitude. 



B. The Ribault phase of the Bechgaard salts 

The organic conductors of the Bechgaard salts family 
(TMTSF)2X (where TMTSF stands for tetramethyhe- 
traselenafulvalene) are well known to have remarkable 
properties in a magnetic field. In three members of this 
family (X=C104, PFe, Re04), a moderate magnetic field 
of a few Tesla destroys the metallic phase and induces 
a series of jSDW phases separated by first-order phase 
transitions .qO 

Accor|ding to the so-called quantized nesting model 
(QNM),eI the formation of the magnetic-field-induced 
spin-density waves (FISDW) results from competition 
between the nesting properties of the Fermi surface and 
the quantization of the electron motion in a magnetic 
field. The formation of a SDW opens a gap, but leaves 
closed pockets of electrons and/or holes in the vicin- 
ity of the Fermi surface. In the absence of a mag- 
netic field, these pockets are too large (due to imper- 
fect nesting) for the SDW phase to be stable. In the 
presence of a magnetic field H, they become quantized 
into Landau levels (more precisely Landau subbands). 
In each FISDW phase, the SDW wave vector is quan- 
tized, Qat = {2kp + NG, Qy) with A'' integer, so that 
an integer number of Landau subbands are filled. [Here 
G — eHb/Ti ^ kp and — e is the electron charge.] As a 
result, the Fermi level lies in a gap between two Landau 
subbands, the SDW phase is stable, and the Hall conduc- 
tivity is quantized :_apy — —2Ne'^/h per one layer of the 
TMTSF moleculesBEl As the magnetic field increases, the 
value of the integer N changes, which leads to a cascade 
of FISDW transitions. The QNM predicts the integer N 
to have always the same sign. While most of the Hall 
plateaus are indeed of the same sign, referred to as pos- 
itive as convention, a negative QHE is also observed at 
certain pressures (the so-called Ribault anomaly) Bu The 
most commonly observed negative phases correspond to 
N = -2andN = -4. 

In the Bechgaard salts, a weak dimerization along the 



chains leads to a half-filled band. Umklapp processes 
transferring Akp = 27r/a are allowed. Thus the forma- 
tion of a SDW at wave vector Qn = {2kp + NG,Qy) 
will be accompanied by a second SDW at wave vector 



Q 



N 



{Akp,0)-Q 



N 



{2kf 



N- 



Vyj, 



there 



is coexistence of phases N and — iV.Ell Note that in our 
notation, Qjy has both the signs of N and Qy reversed 
compared to Qjv- As discussed in the preceding section, 
actual coexistence may occur only if Q^ (like Qn) is a 
good nesting vector. This is the case with the Fermi sur- 
face shown in Fig. Pa (since Qy ^ 7r/&), but not with the 
one shown in Fig. ^a (since Qy ~ -k /2h). 

In Ref. ^ we have studied the effect of umklapp scat- 
tering on the FISDWphases, starting from the Fermi 
surface shown in Fig. |l|a. We have shown that for weak 
umklapp scattering Qy ^ ir/b. In that case, the SDW 
with negative quantum number ~\N\ has a vanishing 
amplitude and can be ignored. However, for N even, 
there exists a critical value of the umklapp scattering 
strength above which the system prefers to form two 
transversely commepsurate SDW's {Qy = i^/b). For cer- 
tain dispersion lawj^ the SDW with negative quantum 
number — |A^| has the largest amplitude, which leads to 
a negative Hall plateau (Fig. ||a). Since the umklapp 
scattering strength is sensitive to pressure, we have sug- 
gested that this provides a natural explanation for the 
negative QHE (Ribault anomaly) observed in the Bech- 
gaard salts.E^ [In the following, the 'negative' phases are 
referred to as 'Ribault' phases.] 

It should be noted that this explanation relies on a 
simple Fermi surface (Fig. ||a) , which does not necessar- 
ily provide a good approximation to the actual Fermi 
surface of the Bechgaard salts. With the more realis- 
tic (according to band calculations) Fermi surface shown 
in Fig. IpD, umklapp processes have, only a small effect 
and do not lead to negative QHEo Therefore, our ex- 
planation of the Ribault anomaly should be taken with 
caution. For it to be correct, the parameter k should be 
smaller (typically not larger than tt/IO) than the value 
7r/4 predicted by band calculation. In the following, we 
consider only the case k = 0. 

We have shown in Ref. |l] that the SDW's in the Rib- 
ault phase are likely to become circularly polarized (he- 
licoidal SDW's) when the umklapp scattering strength is 
further increased (Fig. ||b). The QHE vanishes in the he- 
licoidal phase. We will see that the circular polarization 
also affects the collective modes. 



C. Outline of the paper 

In the next section, we introduce the effective Hamilto- 
nian describing the FISDW phases. The partition func- 
tion is written as a functional integral over a bosonic aux- 

we 



iliary field that describes spin fluctuations. In Sec. Ill 



field propagators and the mean-field particle-hole suscep- 
tibilities. In Sec. [V, we derive the low-energy effective 
action of the SDW phase by taking into account fluctua- 
tions of the bosonic auxiliary held around its saddle-point 
value. We consider only 'phase' fluctuations, i.e. sliding 
and spin-wave collective modes. We do not study ampli- 
tude collective modes, which are gapped and do not cou- 
ple to phase fluctuations in the long- wavelength limit .El 

We find two sliding modes: a (gapless) Goldstone mode 
corresponding to a sliding of the two SDW's in opposite 
directions (out-of-phase oscillations), and a gapped mode 
corresponding to in-phase oscillations (Sec. M). The real 
part of the conductivity exhibit two peaks, which reflects 
the presence of two sliding modes. The low-energy mode 
carries only a fraction of the total spectral weight Loi/A 
(ujp is the plasma frequency) , which is determined by the 
ratio of the amplitudes of the two SDW's (Sec. |V]| ) . 

The spin- wave modes are studied in Sec. VII . There is 



a Goldstone mode corresponding to in-phase oscillations 
of the two SDW's, and a gapped mode corresponding to 
out-of-phase oscillations. The spectral function Imx™* 
is computed in Sec. VIII (x"^*^* is the retarded transverse 
spin-spin correlation function). Both gapped modes are 
found to lie above the mean-field gap in the case of the 
Bechgaard salts. 

In Sec. IX, we study the collective modes in the heli- 
coidal phase where the two SDW's are circularly polar- 
ized. For a helicoidal structure, one cannot distinguish 
between a uniform spin rotation and a global transla- 
tion, so that there are only two phase modes. The Gold- 
stone mode does not contribute to the conductivity and 
is therefore a pure spin- wave mode. Thus all the spectral 
weight in the conductivity a{ui) is pushed up above the 
mean-field gap. Both modes contribute to the spin-spin 
correlation function. 

It should be noted that the long-wavelength modes are 
not the only modes of interest in the FISDW phases. 
There also existj—magneto-rotons at finite wave vectors 
{qx = G, 2G, ...).EaO These modes are not considered in 
this paper. 

In the following, we take fi — ks — ^. 



II. MODEL AND EFFECTIVE HAMILTONIAN 

In the vicinity of the Fermi energy, the electron disper- 
sion law in the Bechgaard salts is approximated as 



E[kx, ky) = VF{\kx\ - kp) + t±{kyb), 



(2.1) 



where k^ and ky are the electron momenta along 
and across the one-dimensional chains of TMTSF. In 



perform a saddle-point approximation, thus recovering 
the mean-field results of Ref. W]. We obtain the mean- 



Eq. (2.1), the longitudinal electron dispersion is lin- 
earized in kx in the vicinity of the two one-dimensional 
Fermi points ±kp, and vp = 2ataSui{kpa) is the corre- 
sponding Fermi velocity (i^ being the transfer integral 
and a the lattice spacing along the chains). Whenever 
necessary, we will impose a ultraviolet energy cutoff Eq to 
simulate a finite bandwidth. The function t±(u), which 



describes the propagation in the transverse direction, is 
periodic: t±{u) = t±{u + 27r). It can be expanded in 
Fourier series as 



t±{u) — — 2ibCos(ii) — 2t26Cos(2M) 

— 2i3b cos(3w) — 2i4b cos(4m) ■ 



(2.2) 



If we retain only the first harmonic {ti,), we obtain a 
Fermi surface with perfect nesting at (2fci?,7r/6). The 
other harmonics t2b, t^b ■ ■ ■ "^ib generate deviations from 
perfect nesting. They have been introduced in order to 
keep a reahstic description of the Fermi surface despite 
the hnearization around rkkp. In the foUowing, we shaU 
retain only tb, ^26 and i4f,. i3b|-does not play an impor- 
tant role and can be discarded.tll We do not consider the 
electron dispersion in the z direction, because it is not 
important in the following (its main effect is to introduce 
a 3D threshold field below which the FISDW cascade is 
suppressedO) . 

The effect of the magnetic field H along the z direc- 
tion is taken into account via the Peierls substitution 
k -^ — iV — eA. (The charge e is positive since the ac- 
tual carriers are holes.) In the gauge A = (0, i/x,0) we 
obtain the non-interacting Hamiltonian 

Ho = ^ / d^r-0t^(r)[wi^(-iaa^ - kp) 

a, a "^ 

+t^{-ibdy-Gx)+ah]^ija,^{r). (2.3) 

Here ipaa-iY) is a fermionic operator for right (a = -I-) and 
left (a = — ) moving particles, a = + (— ) for up (down) 
spin. We use the notation r = {x, mh) {m integer) and 
J d^r = b^^J dx. G = eHb is a magnetic wave vec- 
tor and h = ^bH is the Zeeman energy (we assume the 
electron gyromagnetic fact or t o be equal to two). Diago- 
nalizing the Hamiltonian (2^) we obtain the eigenstates 
and eigenenergies 



0£(r) 



1 



J\!.-T+i{a/i^^)TA_(kyb-Gx) 



yjLxLy 

eQa(k) = €aa{kx) = vpiak^ - kp) + (Th, (2.4) 

where L^Ly is the area of the system, Wc = vpG, and 



T^{u) 



du't±{u'). 



(2.5) 



In the chosen gauge the energy depends only on k^, i.e. 
the dispersion law is one-dimensional. This reflects the 
localization of the electron motion in the transverse di- 
rection, which is at the origin of the QNM (see Ref. pi|for 
a further discussion). Note that contrary to ky, the quan- 
tum number kx is not the momentum since the operator 
dx does not commute with the Hamiltonian. 

The interacting part of the Hamiltonian contains two 
terms corresponding to forward (52) and umklapp ((73) 
scattering: 



«w = f E /<''''*'L(r)*L(r)V's,.(r)*„(r) 



-11: ^ 



ff3v- / ^2^e— 4fe--^t^(r)v;L(r)^„^(r)V'a<.(r), 



(2.6) 



where a = —a and a = —a. For repulsive interaction 
92,93 > 0. We shall assume that umklapp scattering is 
'weak': 173 < 92. We do not consider backward scatter- 
ing (gi) since it does not play an important role in the 
following. The interaction strength is best parametrized 
by dimensionless coupling constants 52 = 92N(0) and 
53 = 53-^(0), where A^(0) = l/nvpb is the density of 
states per spin. 

In a mean-field theory, the cutoff Eq is of the order 
of the bandwidth ta. It is however well known that 
mean-field theory completely neglects fluctuations and 
cannot be directly applied to quasi-lD systems where 
the physics, at least at high temperature, is expected to 
be one-dimensional. The general wisdom is that there is 
at low temperature a-dimensional crossover from a ID to 
a 2D (or 3D) regime.cZlc3 In Bechgaard salts (at ambient 
pressure), this dimensional crossover occurs via coherent 
single-particle interchain hopping. Although the exper- 
imental assignment of the crossover temperature is still 



under debate (see for instance Refs. 29|~pl[), the success 
of the QNM in explaining the phase diagram of the com- 
pounds (TMTSF)2PF6 and (TMTSF)2C104 provides a 
clear evidence of the relevance of the Fermi surface at 
low temperature. In the 3D regime, a mean-field theory 
is justified provided that the parameters of the theory are 
understood as-effective parameters (renormalized by ID 
fluctuations). Ell Eq is then smaller than the bare band- 
width (~ ta) and corresponds to the renormalized trans- 
verse bandwidth. For the same reason, the interaction 
amplitudes 172 and 33 can be different than their bare 
values. The Hamiltonian [Eqs. ( |2.3D and (2.6)] should 
therefore be understood as an effective low-energy Hamil- 
tonian. 

Since collective modes are best studied within a func- 
tional integral formalism, we write the partition function 
as ^ = J 'D^*'D'il} e~^°~^'°* where V^*^ is an anticom- 
muting Grassmann variable. The actions 5*0 and S'int are 
given by 

^0= [d^{j2 /rf'''C.(r,T)9,V^a.(r,T)+7io[^*>]}, 

Sint^-l Y. j d'r dTO*^„{v,T)gc.^,[v)0^,„{v,T), (2.7) 

a^Oi' .a 

where the imaginary time r varies between and /3 = 
l/T. We have introduced the composite field 



Oa^(r, r) = V'Sff(i", T)V'a<T(r, r) 
and the matrix 



(2., 



5(r) 



52 
536 



-i4kFX 



536 



lAkpx 



52 



(2.9) 



Introducing a complex auxiliary field AQo-(r, t), we de- 
couple Si-at by means of a Hubbard-Stratonovitch trans- 
formation. This leads to the action 

S = Sa+ f d^rdTA\{r,T)g{r)A^{r,T) 

]Jd\dTAUr,r)g{r)0,{r,T). (2.10) 



E 



We use the spinor notation Ao- — (A+o-, A_cr)"^ and 
Oc = (0+(j,0_o-)"^. Since Aqc, couples to the 
field OQ(j(r,T) = OQ^(r,T), it satisfies th e con straint 
Aacr(i", t) = Agg.(r, r). Note that the action ( 2.10 ) main- 
tains spin-rotation invariance around the magnetic field 
axis. In the FISDW phase, because of the Zeeman term, 
the magnetization is perpendicular to the magnetic field 
axis, so that the only spin-wave mode corresponds to 
rotation around the z axis. In zero magnetic field, a dif- 
ferent approach should be used in order to maintain the 
full SU(2) spin symmetry (see Ref. 0). 



III. MEAN-FIELD THEORY 

In this section we look for a mean-field solution cor- 
responding to a phase with two sinusoidal {i.e. linearly 
polarized) SDW's: 



{Sx{r)) = ^ MpN cos{(j)pN) cosiQpN ■ r + 9. 



pN) 



p=± 



{Sy{r)) = ^ MpN sin(0pAr) cos(QpAr • r + Op^) 

p=± 
{S.{r)) - 0, (3.1) 

where S^{r) = 'Eaaa' 'l^*aA'^Ka'^c,cT'{r) is the spin- 
density operator and t'^ {v = x, y, z) the Pauli matri- 
ces. Because of the Zeeman coupling with the magnetic 



J 



field, the SDW's are polarized in the [x, y)-plane. The 
variable (j)p]si determines the polarization axis, while OpN 
gives the position of the SDW's with respect to the un- 
derlying crystal lattice. 

We assume that the external parameters (magnetic 
field, pressure...) are such that the system is in the Rib- 
ault phase, characterized by a negative QHE and the 
coexistence of two SDW's with comparable amplitudes. 
We choose the sign of N such that A'' refers to the SDW 
with the largest amplitude [AIm > ^n)- [^ is even and 
negative in the Ribault phase.] 

The mean-field solution corresponds to a^addle-point 
approximation with a static auxiliary fielcS 



A„,(r) = (0„,(r,r)} 



^^(pjV)gmQp„-r^ (3.2) 

p 



The relation between Aaa and {Oaa) results from the 
stationarity condition of the saddle-point action (see 
Eq. (3.4) below). Because of the constraint A*g.(r) = 

Aag^(r), the order parameters satisfy Aaa = Ag^ ' . 
Among the eight complex order parameters, only four 
are therefore independent and sufficient to character- 
ize the SDW phase. The order parameters AcTo- = 
Aqct le"^""" are related to the spin-density modulation 



]y)by 



|AiP^)| = 



M, 



pN 






'A.^=-^(^?r+v'LtO- (3-3) 



Since by convention Afjv > Mf^, we have |A 



\Aaa\. Note that the condition Aao- '* 



Wi 



ipN) 






> 



Ag^ imposes 



The mean-field action is given by 



5MF=/3y'dVAJ(r).9(r)At(r)-^y'dVdTy'd2r'dr'(CT(r,r), VSi(r,T)) 

(5(r-r')<5(r-r')A„t(r)\ fi^^^ir^r') 



i(0)-l^ 



G':r\r,T;r',T') 



5(r-r')<5(r-r')A;^(r) 



-.(0)/ 



Here G)^J{r,T;r' ,t') is the Green's function in the absence of interaction (52 = 53 
AQ,cr(r) acting on the electrons is defined by 

A„,(r) = 52Ag,(r) + g3e^"^'""Aa<,(r), 



(3.4) 
0). The effective potential 

(3.5) 



i.e. A(j(r) — g{r)Acr{r). Rewriting the action in the basis of the eigenstates ^g of Wo [Eq. ( |2.4|) ], we obtain 
SMF = (3jd\A\{r)g{r)A^{r) 



a.k.k' 



(V':^(k,c.), Ci(k,c.)) 



Aai(k,k') 



A„T(k,k') 
'5k,k'(«w-egj,(fc2;)) 



■0a|(k',w) 

■0ai.(k',w) 



(3.6) 



We denote by a; = 27rT(n + 1/2) (n integer) fermionic 
Matsubara frequencies. We emphasize that here k and 
k' refer to the quantum numbers of the eigenstates (f>^ 
of Hf). We have introduced the (particle- hole) pairing 
amplitudes 



A„,(k,k') = / d'rC*(r)C'(r)A„,(r) 



n— — CO 



J ^ian{kyb-n /2) ^ 



(3.7) 

where /„ = In{% — ^/^)- The coefficients In{qy), defined 
by 



Iniqy) 



2ir 



^ inu+{i/LJa)[T±{u+qyb/2)+T^{u~qyb/2)] 



2lT 



(3.8) 



are well know in the QNM. They depend on the trans- 
verse dispersion law t^ikyb) and measure the degree 
of perfect nesting of the Fermi surface. For a per- 
fect nesting at wave vector (2kp,Tr/b), one would have 
In{T^/b) — Snfi, while JH the generic situation (imperfect 
nesting) |J„(gy)| < lS3 The convention Mj^ < Mn im- 
plies |/jv| < l-^Afl- Note that k^ is coupled not only to 



k^ - aQT"' but also to k^ - aQT"' + anG [Eq. (|37D]. 
This reflects the fact that k^ is not the momentum in the 
gauge A = {0,Hx,0). 

From now on, we shall use the quantum limit ap- 
proximation (QLA) (also known as the single gap 



approximatioiu). This approximation, which is valid 
when LUc is much larger than the temperature and the 
SDW gap, consists in retaining only the most singu- 
lar (particle-hole) pairing channels. In the metallic 
phase, these channels present the logarithmic singular- 
ity ~ \n{Eo/T) (which is characteristic of a SDW system 
with perfect nesting). This singularity results from pair- 
ing between electron and hole states of the same energy 
and is responsible for the opening of a gap at the Fermi 
level. While the QLA strongly underestimates the transi- 
tion temperature and the value of the gap at low temper- 
ature, and also neglects gaps opening above and below 
the Fermi level, it is an excellent approximation when 
calculating -the collective modes (we shall come back to 
this point ).E3 

Thus, in the QLA, Accr(k, k') is not zero only if 
^aa{kx) = —('as{k'x), which requircs k'^ ^ k^ — a2kp. 
Only the term with n = pN survives in Eq. (3.7), which 
yields 



'Qcr(,iijK j — dk',k— qQo^qctV"'!/ ji 

A„,(fc,)=.^A(r)e-''^('=«^-'^/2), 



(3.9) 



where Qo = {2kp,7T/b). We have introduced the new 
order parameters Aaa related to Aaa by 



A(rJ^) = /p^(g2A(r)+53Ar^) = A 



ipN)* 



A(p^) = 



g2lpNA.aa - g^IpN Ag^„ 

{92-93)InIn 



(3.10) 



J 



In the QLA, the mean-field action therefore reduces to a 
2x2 matrix (as in BCS theory): 



Smf = (3 I d\A\{r)g{r)AT(r) 



E 

a,k,c 



iuj-ea-iikx) Aa]{ky) 

A*^^{ky) iuj + ea](kx) 



■0Q|(k,u;) 
'0ai(k-aQo,tj) 



(3.11) 



using ea[{kx — a2kp) = —ea'\{kx)- It is remarkable 
that within the QLA the mean-field action can still be 
written as a 2 x 2 matrix. The presence of a second SDW 
changes the expression of the pairing amplitude (which 
becomes fcy-dependent), but not the fact that the state 
(k, — ) couples only to (k -I- Qo,-f). This is not true in 
zero magnetic field where the presence of a_aecond SDW 
leads to complicated mean-field equations.Ej 



r 



En, where T 
state energy: 



and En — ~N{0)Eq is the normal 



AE ^ 









ipN 



iV(0) 



E 

p 



A 



{pN),2 
aT 



^(0)|X(iV)|2 
2 I'^aT I 




(3.12) 



A. Ground-state energy 

The mean-field action being Gaussian, we can integrate 
out the fermion fields to obtain the ground-state conden- 
sation energy (per unit area), AE = —{T\nZ)/LxLy — 



To obtain ( 3.12 ), we have explicitly taken into account 

the fact that lAi^^^l > lAi^^^l and \AiV\ > \AiV\. The 
latter inequality follows fr om M ^ > Mff, while the for- 
mer is derived below Eq. (3.15). Using the gap equation 

dAE/dA^^ = 0, we can rewrite the condensation en- 
ergy as 



AE=- 



NiO) 



E 



I a\ I 



(3.13) 



a form which is reminiscent of the BCS-Ukc expression of 
a SDW system with perfect nesting (A£; = -7V(0)|A|2/2 
with A the SDW gap). As stated above, the QLA 
amounts to taking into account only the 'perfect nest- 
ing component' of the particle-hole pairing. 

The requirement that AE is minimal with respect to 
variation of the order parameters leads to constraints on 



the phases of the complex order parameters Aq 



(pN) 












if InIn > O7 
TT if InIn < 0. 



(3.14) 



For tib = 0, /^ - (-l)^/2/Ar. 

smah, sgn(/Ar/^) = (-1)^/^. 
(3.14), we then obtain 



When tib is finit e but 
Using Eqs. (|3.10|) and 



AN) 



|a(4Va^T^| 



AiNJi. 



N\ 



A-- |7^//^|-r|A^-VA 



(3.15) 



qT 



Eq. ( p5| ), together with \A^^^ j A^^„\ = A%/M^ < 1 



W/a(^)i 



and |/^//jv| < 1, implies IAq^YAqo- | < 1. It also shows 



.W/aW 



(JV) 



that Ag^'/Ajjt is a real number, so that the phases (^^^ 

'at 



and ip-J are identical modulo tt: 



(Af) (TV) r 1 



From (3_^) and ( 3.16| ), we then deduce 



<iN 



^Af; 



TN = 



'N- 



(3.16) 



(3.17) 



In the mean-field ground state, the two SDW's have 
the same polarization axis {(pN ~ 4'n)- The condition 
Ofq = —0f[ means that the two SDW's can be displaced 
in opposite directions without changing the energy of the 
system. They can therefore be centered either on the lat- 
tice sites or on the bonds. For 9n = Off = 0, both SDW's 
are centered on the sites, while for 9n = ^(^n — '"'/^j they 
are centered on the bonds. The condition 9^ + 0^ = 
is related to the pinning that would occur for a commen- 
surate SDW. Indeed, for a single SDW with wave vector 
{2kp, K/b), it becomes 29 = where 9 is the phase of the 

sDw.y 

This degeneracy of the ground state results frorpj'ota- 
tional invariance around the z-axis in spin space ,0 and 
translational invariance in real space. The latter holds 
in the FISDW states, since the SDW's axe incommensu- 
rate with respect to the crystal latticcEj According to 
the mean- field analysis, and in agreement with general 
symmetry considerations, we therefore expect two (gap- 
less) Goldstone modes: a spin-wave mode corresponding 
to a uniform rotation around the z-axis of the common 



polarization axis, and a sliding mode corresponding to a 
displacement of the two SDW's in opposite directions. 

Without loss of generality, we can choose the phases 
(faa and ifaa cqual to or n. The order parame- 
ters Aaa and Aaa are then real (positive or negative) 
numbers. Furthermore, since the constraints obtained 
by minimizing the free energy relate (fan (}Paa- ) only 



AN) 



to <^acr (<Pgo- )j we can also impose (pa 
(Paa . With such a choice. 



(N) 






Lp\^J and 



(3.18) 



are real (positive or negative) numbers independent of a 
and a. 9pN = and (J)n = 4'n = or tt. The SDW's are 
polarized along the x axis. Introducing the notations 



53 



52 



, c 



'AT 



A 



7 = 



N 



In' 

_An 
An' ' An' 

(with I7I = Mj^/Mn), we deduce from (|3.10|) 

7 -|- r 7 ~ '"C 



7 = C 



1 -|- r7 ' 



7 



C-n' 



(3.19) 



(3.20) 



Note that 7 and C have the same sign [see ( p.ld^ , and 
I7I, \l\,r, Id < 1. In the Ribault phase, I7I ~ ItJEI Both 
SDW's have the same amplitude when |C| = 1- In that 
case 7 = C and 7 = 1. 



B. Mean- field propagators 

The normal and anomalous Green's functions are de- 
fined by 

Gaair,r',uj) ^ -{iljaair,uj)ijl„ir',uj)) 

= ^C(k,k',c.)0£(r)C*(r'), 

k,k' 

= ^i^„.(k,k',c.)C(r)C* (!■')• (3.21) 

k,k' 



From (3.11), we deduce that 

G'Q<x(k, k',w) = (5k,k'GQo-(k,w), 

i^QCT(k, k',w) =(5k',k-aQo^a<T(k,w), 

with 



(3.22) 



GaaiK^) 



FaaiKu}) 



'iuj-eaa{kx) 






(3.23) 



The excitation energies are given by ±_Bcto-(k), where 



Writing Aif ) = |A„^ 



;L(fc.) + |A„.(fc,) 



(3.24) 
we obtain (for N even) 



p 
X cos(2a7Vfc,6 + ^(fj - <^W). (3.25) 

The quasi-particle excitation gap 2\Aaa{ky)\ depends on 
the transverse momentum ky (Fig. 0). When both SDW's 
have the same amphtude, IAatJ = |Ajy|, the spectrmxi be- 
comes gapless since A^aiky) vanishes at some points on 
the Fermi surface. 

We have shown in Ref. ^ that the occurrence of heh- 
coidal phases prevents the spectrum from becoming gap- 
less. Indeed, the stabihty of the sinusoidal phase requires 



mmkJAaaiky)\ 



1 - lo-l 

- — P^ > 0.32, 

l + l7l 



(3.26) 



maxkJAaa{ky)\ 

i.e. I7I < 0.52. Eq. ( ^.26 ) is deduced from a Ginzburg- 
Landau expansion of the free energy (valid only close to 
the transition line) and may be slightly changed at low 
temperature. 

In the helicoidal phase, there is no dispersion of the gap 
with respect to ky (see Sec. px|) . Since the spectrum be- 
comes gapless for I7I = 1, it is natural that the system. 



J 



above a certain value of I7I, prefers to form helicoidal 
SDW's in order to lower the free energy by opening a 
large gap on the whole Fermi surface. 



C. Mean- field susceptibilities 

The results reported in this section are derived in detail 
in Appendix IaI We study the mean-field susceptibilities 
(which will be useful in Sec. |rv| ) 

Xaa,a'a'ir,T;r',T') = {Oaair , t)0*^, „, {r' , t')) 

-(0„.(r,T))(0;,,,(r',r')), (3.27) 

where the mean value (• • •) is taken with the mean- field 
action 5mf [Eq. ( |3.1l| )]. Since the latter is Gaussian, we 
immediately deduce the only two non-vanishing compo- 
nents 






'Gaa{r, r; r', r')Gaff(r', r'; r, r), 

-^„,(r,r;r',r')F„,(r',r';r,^)- 

(3.28) 



At T = and for qy = 0, the Fourier transformed 
susceptibilities are given by (Appendix ^) 



x+T.+T(q + Qpw= q + QpN, (^u) 
x+T -i(q + QpN, q - QpN, uju) 

X+T,-i(q + QpN, q - QpAT, ujv) 



iV(0) 



m^L 



^N\ 



■vlql 



6(1-72)A 



InInI 



^Wx 



,,.,ffl)4(i_,^, 



3(1-72)A 



■'P.+- 



N{Q) 



1 



vlql 



6(1-7^)A^ 



InInI, 



(3.29) 



where uji 
frequency. 
v1,ql,ujl < 



= 2itTv 
Eqs. (I 

A^(l 



V integer) is a bosonic Matsubara 
) are valid in the low-energy limit 
7^). We have used the QLA (see 
Appendix |AJ) , whose validity is discussed in detail in Ref. 
p5| . Although the QLA does not predict accurately the 
value of the mean- field gaps Ap^ and Ap^r, it is an ex- 
cellent approximation of the low-energy properties (such 
as the long-wavelength collective modes) once the values 
of ApN and ApM are known. Since the absolute value of 
the latter do not play an important role for our purpose, 
the QLA turns out to be a perfect mean to compute the 
collective modes. 



IV. LOW-ENERGY EFFECTIVE ACTION 

In this section we derive the low-energy effective ac- 
tion determining the sliding and spin- wave modes. We do 
not consider amplitude modes which are gapped and do 
not couple to 'phase' fluctuations in the long-wavelength 
limit. We consider only longitudinal (i.e. with the wave 
vector parallel to the chains) fluctuations {qy — 0). These 
fluctuations are of particular interest since they couple 
to an electric fleld applied along the chains (which is a 
common experimental situation): see Sec. VI on the op- 
tical conductivity. Including a flnite qy would allow to 
obtain the effective mass of the collective modes in the 
transverse direction. Such a calculation is however much 
more involved and will not be attempted here. 

The low-energy effective action is derived by study- 



ing fluctuations of the auxiliary fleld Aco-lr, t) around 
its saddle-point value Aaa{j^)- We write AQ,cr(r, r) = 
Aqct (r) + rjaa (r, r) and calculate the effective action to 
quadrat ic or der in the fluctuating field rj. The fermionic 
action (2.10) can be rewritten as 



^[V'*>] = 5mf[^*>] -E/'^'^'^^^i(^'^)^-('"'^) 



+ / cfrdT[r]\{r,T)g{r)T]^{r,T) 
+ (r/|(r,T)At(r)+c.c.)]. 



(4.1) 



Integrating out the fermions, we obtain to quadratic or- 
der in 77 the effective action 



S[v 






a.a' .a.a' 



Sa^a'Si'C - v')6{t - T')gaa' (r) 

XcKT,QV'(r,T;r',r') r]a'a'{r' ,t'), (4.2) 



where the susceptibility x is defined by ( A13| ). There is 
no linear term since we expand around the saddle point. 
Considering phase fluctuations only, 



A„<,(r,T)=^Apwe^"'^- 



iV-r+ivL",, '(r,r) 



(4.3) 



(ySg^ (r, t) = — ((3qct (r, t). Using this latter condition, 
we deduce from (|4.2|) the effective action for the phase 



(4.5) 



X \g2Sa,a'Sp.p' + g?i5a,a'5p,p' + ^^q/ (?) : 

where q— [q^jUJi,), and 



(4.6) 



In Eq. ( |4.6| ) q = (g^;, 9y = 0). Using Eq. ( [A15| ), and in- 
troducing charge and spin variables as in the mean-field 
theory [see Eq. (3.3)], 



0,^.(r,r) = i(^V'r(r,r)-^Lt)(r,r)), 

</)p^(r,r) = -i(^V'r(r,r)+^Lf^(r,r)), (4.7) 

we obtain a decoupling of the sliding and spin-wave 
modes, which are determined by the actions 



1 



S[0] ^-J2ie^{~q), 0^{-q))V-^{q) 



riq) 



^ ) ' (^-^^ 



we have 

,7„. (r, r) = * ^ A.^e-Q''" ^^i^f^ (r, r) 



(4.4) 



to lowest order in the (real) phase field ifaa (r, r) . 
Note that the property A*g.(r, r) = Aaa{^,T) imposes 



J 



^c"h(9)=492A 



S[<P] -lj2i^N{-q), cbj,i^q))V-^\q) ( J^[| 



1 — C-I-+ — r^c + 2rc-\ 

^ ^ -7r(l - C++ - c ) - 7(1 -I- r^)c^ 



^7p^(^1=452A^( ;,(i't 



/ 1 — C++ — r^c — 2rc4 



C++-C— )- 7(1 + H)C4 



The coefficients Cpp' = Cppi{q) are linear combinations 
of the mean- field susceptibilities x ^^d are defined in Ap- 
pendix |A|. 



with 

-7r(l - C++ - c — ) — 7(1 + r'^)c+- 
7^(1 — c - r^c++) + 27^rc+_ 

7r(l - C++ - c — ) - 7(1 + r^)c+- 
7^(1 — c__ — r'^c+j^) — 27^rc+- 

V. SLIDING MODES 



(4.9) 

(4.10) 
(4.11) 



The effective actions [Eqs. (4.8) and (4.9)] have been 
obtained by expanding the auxiliary field with respect to 
the phase fluctuations ipaa (r, t) [Eq. (4.4)]. While this 



The dispersion of th e slid ing modes is obtained from 
detX»7jj^((?) = 0. From (|t|), we deduce 

AeiV-2{q) = W2^%l\l - r^m - c++)(1 - c-) 



procedure turns out to be correct at zero temperature, 
an alternative and more rigorous approach consists in 
performing a chiral rotation of the fermion fields, which 
allows to directly expand with respect to dr^Paa (r,T) 
and Viy9acr (r,T).aCJ The drawback of this approach is 
that one has to calculate the non-trivial Jacobian of the 
chiral transformation (the so-called chiral anomaly). 



-r'^c++c + 2rc+_ - (1 - r'^)c^_^ 



(5.1) 



The gap equation (A12) shows that dei'D^2{q = 0) = 0, 



which indicates that there is a gapless mode with van- 
ishin g en ergy in the limit of long- wavelength {q^ — > 0). 
Eq. (5.1) simplifies into 



'5C++ — 

7 



Sc —^ — 2r6cj 



c 



-(1 — r"')((5c_|_+(5c_ 



5cl_ 



(5.2) 



A. Goldstone mode 



Eq. ( ^.2| ) admits the solution ul + v^q^ = 0. After an- 
alytical continuation to real ft'cquencies (iui, — > uj + iO^), 
we obtain a mode with a linear dispersion law (Fig. ]^ 






(5.3) 
(5.4) 



As expected from the mean- field analysis (Sec. til), the 
Goldstone mode corresponds to an out-of-phase oscilla- 
tion of the two SDW's. We do not find any renormaliza- 
tion of the mode velocity because we have not taken into 
account long- wavelength charge fluctuations. The latter 
couple to the sliding modes and renormalize the velocity 
to vp{l ~^g4^N{0)y/^, where g^ is the forward scattering 
strength.EJ In order to obtain this velocity renormaliza- 
tion, it would have been necessary to include forward 
scattering in the interaction Hamiltonian and to intro- 
duce an auxiliary field for the long-wavelength charge 
fluctuations.D 



B. Gapped mode 



The second mode obtained from ( |5.2| ) corresponds to 
a gapped mode with dispersion law (Fig. 0) 



2 2 
VfQx 

12 



3 + 57 



o,2 



52 1 



37 



^,2 



AatAjv 



In In 



(5.5) 
(5.6) 



The oscillations of this mode satisfy 
On rC - 7 ^(3 + 7^ + 47C)^ + (1 " ^)^4C7(3 + 7^) 



^N 



n-C 



{3 + j^ + Ar^Cy 



(5.7) 



If we consider the case where the two SDW's have the 
same amplitude (which is not physical because of the 
appearance of the helicoidal phase), which occurs when 
7 = C = ±1, we obtain the simple result 9n — Oj^. 

Using the physical parameters of the Bechgaard salts, 
we find that luq is generally larger than the order pa- 
rameters |A-|-jv|, so that the gapped sliding mode ap- 
pears above the quasi-particle excitation gap (in gen- 
eral within the first Landau subband above the Fermi 
level). We therefore expect this mode to be strongly 
damped due to the coupli ng w ith the quasi-particle ex- 
citations. Note that Eq. (|5.6|) giving wq is correct only 
if ujQ ^ [AatKI — 7^)"^/^, since it was obtained using 



expressions of the susceptibilities Cpp' valid in the limit 



,vhl<^K(.^ 



■i') 



It is worth pointing out that the two sliding modes 
bear sepie similaiities with phase modes occurring in two- 
bandtll or bilayeiila superc ondu ctors. [This also holds for 
the spin- wave modes (Sec. VII) and the collective modes 
of the helicoidal phase (Sec. IX).] 53/52 plays the same 



role as the ratio between the interband (or interlayer) and 
intraband (or intralayer) coupling constants. There are 
also analogies with phase modes in d -I- id' superconduc- 
tors, and, to some extent, with,pJasmon modes occurring 
in conducting bilayer systems.EJ While the correspond- 
ing phase modes in superconducting systems have not yet 
been observed, plasmon modes in semiconductor double- 
well structure have been observed recently via inelastic- 
light-scattering experiments O 



VI. OPTICAL CONDUCTIVITY 

At T = there are no thermally excited carriers above 
the SDW gap. The response to a low-frequency electric 
field is then entirely due to fluctuations of the SDW con- 
densate. Since the electric field does not couple to the 
amplitude modes, we only have to consider the phase 
fiuctuations studied in the preceding sections. In the 
next section, we determine the charge (puw) and cur- 
rent (jDw) fluctuations induced by phase fluctuations of 
the condensate. In Sec. VI B| , we calculate the current- 
current correlation function (jdwJdw) which determines 
the conductivity. 



A. Charge and current operators 

In this section we calculate explicitly the charge in- 
duced by phase fluctuations of the condensate. The cur- 
rent is then obtained from the continuity equation. 

We consider a source field Ao{r,T) which couples to 
the (long-wavelength) charge density fluctuations. This 
adds to the fermionic action the contribution 

A5[^*,^] = /'dVdrAo(r,T)^Ca(r,T)^a.(r,T). 

(6.1) 

The charge operator is then obtained from the action 
by functional differentiation with respect to the external 
field Aq: 



p(r,T-) 



SS 



SAoir,- 



(6.2) 



Ao=0 



Following Sec. ||, we introduce the auxiliary field 
AQ(j(r, t) and integrate out the fermions. This leads to 
the action 



10 



S[A*,A,Ao] = / d\dTA\{r,T)g{r)A^{r,T) 



^Trln(-e-i+io), 



(6.3) 



where 



J 



G':r\r,T;r',T') 



g-Ur T-r' r')= I ""^ 

V'^(r-r')<5(r-r')A;^(r,r) 



(5(r-r')<5(r-r')AaT(r,T 
.(0) 



Gl°^-i(r,r;r',T') 



(6.4) 



io(r, r; r', r') = 5(r - r')5(r - r')Ao(r, r)i. (6.5) 

We denote by 1 the 2x2 unit matrix. Expanding the 
action with respect to Aq, we obtain 

S[A*,A, Ao] - ^[A*, A] + J2 Tri^aio] + 0{Al), (6.6) 



where S'[A*,A] is the action without the source field. 
Here Tr denotes the trace with respect to time, space 
and matrix indices. From (S.2), we then obtain the fol- 
lowing expression of the charge operator 



Pi^^'^) ^^i^Sai^:Ti^^'^)y 



(6.7) 



where tr denotes the trace with respect to the matrix 
indices only (i.e. TrO = J d^r drtr 0(r,T;r,T) for any 
operator O). 

To calculate p, we write 



A,(r,T) = A,(r)+77,(r,r), 
Va{r,T) =.g(r)77^(r,T). 



Then we have 



^QV- I ' 1 '^ J 



S{r-r')S{T-T') 
??qt(i",t) 



(6.8) 
(6.9) 

(6.10) 



where G^^ is the mean-field propagator. The charge op- 
erator is given by 

p(r,r)=^trgr(r,r;r,r) 

+ J2 tr[erS«Cn (r, r; r, r) + 0{f,'). (6.11) 



The first term in the rhs of (6.11) is the uniform charge 
density po in the mean-field state. The other terms corre- 
spond to charge fluctuations pdw induced by the conden- 
sate fluctuations. To lowest order in phase fluctuations, 
we obtain 



Pdw 



;r,r) = 5:tr[C?rS.C?r](r,r;r,r) 



(6.12) 



r 



Considering only phase fluctuations [Eq. ([4.4[)], we ob- 
tain (Appendix pi) 



PDw(r,r) = — 5a;6'jv(r,T), 

TTO 



(6.13) 



where 

0N{r,T) 



~J^) 



{N), 



'f+^ i^,T) -^pl^'{r,T) 



(6.14) 



~(pAf)/j. ^ ^ g2ApNipaf^ (r, r) + g^ApNifiaf^ (r, t) 



'P 



g2ApN + 93 Ap 



pN 



Note that pDw(r, t) depends only on 9m{y, t), i.e. on the 
largest component of the effective potential Aa^{r,T). 
Using the continuity equation (in imaginary time) 



iedrPuwir^r) + (9a;jDw(r,T) = 0, 



(6.16) 



we obtain the current fluctuations induced by the phase 
fluctuations: 



JDw(r,T) 



Trb 



drON{r,T). 



(6.17) 



Note that the current jdw is parallel to the chains. Con- 
densate fluctuations do not induce current fluctuations 
in the transverse direction. 

This simple result [Eqs. ( 6.13| ) and (6.17)] can be un- 
derstood as follows. For a system with a single SDW, 
the charge fluctuation pDW = dxO/irb is obtained by 
requiring the SDW gap to remain tied to the Fermi 
level. Let us briefly recall this argument. A fluctuat- 
ing SDW potential A„<,(x,t) = A^^e*"2fcf^+^¥'»^(^,^) ~ 
A„^e*"2/cf x+ixa,¥'„„(:E,r) couples the state {k^, j) to (fc^-h 
2kp + 9j,<y9+|(a;,r),t), and (k^,^) to (fc^; -|- 2kF - 
dxf-Tix, t), I). [For simplicity, we neglect the transverse 
direction.] The resulting gap will open at the Fermi level 
only if the system modifles its density in such a way that 



SkF+^ix,T) +SkF-iix,T) = dx(p+-\ix,T), 
5kF+i{x,T) +5kF-]{x,T) = -dx'^--\{x,T). 



(6.18) 



Here we denote by a{kF -f 5kFaa{x,T)) the Fermi wave 
vector on the {aa) branch of the spectrum. This Fermi 
wave vector is time and space dependent because of the 
density fluctuations. Eqs. (6. IS) imply the charge density 
variation 
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Sp{x,T) = —-'^\6kF+a{x,T) +(5fcF-cr(a;,r)j 



2Trb 

= l-d.J{x,r). 
no 



(6.19) 



Let us now go back to the case of interest. We should 
consider the effective potential Aaa{^,T), since it is 
AQ,cr(r, r) (and not Aaa{^, t)) which couples to the elec- 
trons. In pr esence of two SDW's, it is not possible to 
satisfy (6.18) for the two components {N and —N) of 
the effective potential. Nevertheless, it is natural to as- 
sume that the largest gap rem ains t ied to the Fermi level. 
This immediately yie lds E g. ( |6.13 ), and, from the conti- 
nuity equation, Eq. ( 3.17 ). [The latter can also be ob- 
tained by a similar argument without usin g the conti- 
nuity equation.ld] In conclusion, Eqs. ( 3.13 ) and ( 6.17 ) 
express the fact that the largest gap remains tied to the 
Fermi level. 

Finally, we can express the charge and current fluctu- 
ations directly in terms of 9n and 9 



'N- 



PDw(r,T) = 



1 9:r^w(r, r) - r-/d.j:9[^{r,T) 



nb 



1 + r7 



■ / ^ ie dr9N{v,T) ~ ndr9fj{r,T) 

jDw(r,T) = — . (6.20) 



irb 



1 + rj 



B. Current-current correlation function 

In this section we calculate the correlation function 

n(w,.) = (JDw(^^)JDw(-Wi.)), (6.21) 

where Jdw('j-'!^) = JDw(q = 0,[x)y) is t he F ourier trans- 
form of the current operator given by ( p. 17 ) . We have 



Ii{LU,) 



(6.22) 



where V^p {lj^) = (6'pw(w^)6'p'jv(-c^^)) [see Eq. (gj)]. 

The conductivity is defined by cr{uj) — i/{Lo + 
iO~^)W^^{uj) where the retarded correlation function 
n''°'(ti;) is obtained from 11(^^1,) by analytical continu- 
ation to real frequency iuji, —> lo + iO'^ . Using the expres- 
sion of Cpp' (Appendix Kh, we obtain the dissipative part 
of the conductivity 



Re[a(^)] = -f 



,, ,3(1-72) ^^ , ^47 



57 



3 -h 572 



where 



^'ie^vp jb is the plasma frequency. 



Eq. (|6.23 ) satisfies the conductivity sum rule 



dwRe[cr(tj)] = -j-. 



(6.24) 



Thus all the spectral weight is exhausted by the col- 
lective sliding modes. Quasi-particle excitations above 
the mean-field gap do not contribute to the (longitudi- 
nal) optical conductivity, a result well known in SDW 
systems. El Because both modes contribute to the conduc- 
tivity, the low-energy (Goldstone) mode carries only a 
fraction of the total spectral weight. We obtain Dirac 
peaks at iwo because we have neglected the coupling of 
the gapped mode with quasi-particle excitations. Also, 
in a real system (with impurities), the Goldstone mode 
would broaden and appear at a finite frequency (be- 
low the quasi-particle excitation gap) due to pinning by 
impurities. In the clean limit, which is appropriate in 
(TMTSF)2X salts, the presence of impurities does not re- 
store any significant spectral weight to quasi-particle ex- 
citations above the mean-field gap.l Therefore, the frac- 
tion of spectral wei ght c arried by the two modes is cor- 
rectly given by Eq. (6.23). By measuring the optical con- 
ductivity <j{lo\ we can therefore obtain the ratio I7I ~ I7I 
of the amplitudes of the two SDW's. 

We have shown in Ref. 21 that the sinusoidal phase 
becomes unstable against the formation of a helicoidal 
phase when I7I reaches a critical value 7c. By compar- 
ing th e me an-field condensation energies of the sinusoidal 
[Eq. (|1|)] and hehcoidal [Eq. (|9j) in Sec. g below] 
phases, we find 7c ^ 0.4, which is in good agreement with 
the result 7c ^ 0.52 obtained from the Ginzburg-Landau 
expansion of the free energy close to the transition line.Eil 

Since I7I varies between and 0.4 in the sinusoidal 
phase, the Goldstone mode can carry between ~ 100% 
and ~ 50% of the total spectral weight (Fig. ||). When 

this fraction is reduced below ~ 50% (i.e. when I7I ^ 
0.4), a first order transition to the helicoidal phase takes 
place. In the helicoidal phase, there is only one Goldstone 
mode (see Sec. IX). This mode is a pure spin- wave mode 
and does not contribute to the optical spectral weight. 

Experimentally, the helicoidal phase can be stabi- 
lized by decreasing-jpressure, which increases the ratio 
r = 53/52 (Fig. ||).Ej Fig. ^ shows I7I and the low-energy 
spectral weight as a function of r for iJ ~ 10 T. In the 
helicoidal phase, 7 is defined by Eq. (9 



VII. SPIN- WAVE MODES 



For the spin- wave modes, the condition Ae\X>^^{q) ~ 



yields 



(6.23) 



^c++ V 6c —— + 2r5c^ 

7 C 

-(1 - r'^){6c++6c— - 6c\_) = 0. 



(7.1) 
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A. Goldstone mode 



From Eq. (7.1) we deduce the existence of a mode with 
a hnear dispersion law 



w = vpqx, 
l)N[qx,uj) = cj)f![qx,uj). 



(7.2) 
(7.3) 



The spin-wave Goldstone mode corresponds to in-phase 
oscillations of the two SDW's, in agreement with the 
conclusion of the mean-field analysis (Sec. |l|). We do 
not find any renormalization of its velocity because we 
have not taken into account the coupling with the long- 
wavelength spin fluctuations]^ 



B. Gapped mode 



The other solution of (7T) corresponds to a gapped 
mode with the dispersion law 



u:^-- 


2 2 2 


. ,2 


12 r 1-72 


"^1 - 


32 1 - r-^ 7^ 



A«A 



N^N 



In In 
The oscillations of this mode satisfy 

(f>N{qx,i^) = -n(l>^{q^,uj). 



(7.4) 



(7.5) 



In the gapped made, the oscillations of the two SDW's 
are out-of-phase.t3 As expected, the largest SDW has 
the smallest oscillations. Like the gapped sliding mode, 
this mode is found to lie in general above the quasi- 
particle excitation gap, and is therefore expected to be 
strongly damped due to coupling with quasi-particle ex- 
citations. Thus, the spin-wave modes are similar to the 
sliding modes as shown in Fig. y. 



VIII. SPIN-SPIN CORRELATION FUNCTION 

In this section we calculate the transverse spin-spin 
correlation function. For real order parameters Aaa , 
the mean- field magnetization (S(r)) is parallel to the x 
axis. Transverse (to the magnetization) spin fluctuations 
corresponds to fluctuations of the operator 



Syir.T) = -^crO„^(r,T). 



(8.1) 



The correlation function Xyy = (SySy) can be expressed 
as a functional derivative of the free energy with respect 
to an external field that couples to the spin operator Sy. 
This standard procedure allows to rewrite Xyy ^s 



X,,(r,r;r',r') = J ^ aa'(A„.(r, r)A:,,,(r', r')) 

a, a', (7,(7"' 

-i5(r~r')5(r-r')^ff-i(r), (8.2) 

a.a' 

where the mean value (■ • •) is taken with the effective 
action S'[A*, A]. In the following, we drop the last term 
which does not contribute to the spectral function Imx™*. 
To first order in phase fluctuations 

Xyy{r,T;r',T') = 4^ ApArAp'Arcos(QpAr • r) 
p,p' 

xcos{qp,N-r')VPP'{r,T;r',r'), (8.3) 

where 2?PP'(r,r;r',T') = (0pw(r,T)0p,jv(r',T')). Taking 
the Fourier transform, wc obtain 

Xyyi^l + CtQpN, q + a'Qp'AT, CJi.) = ApA? Ap/ArX'P^ (q), 

(8.4) 

where q = {qx,qy = 0). Because of the presence of 
SDW's, the spin-spin correlation function is not diagonal 
in momentum space. In Eq. (3.4), q corresponds to the 
momentum and energy of the spin-wave mode and tends 
to zero for long-wavelength fluctuations. We therefore 
consider the spectral function ImTrqX™*, where Trq is 
a partial trace corresponding to a given spin-wave mode 
momentum q ~ [q^ ,qy — 0): 

TrqXyy = ^ Xyy{(l + aQpN, q + aQpN, i^v) 

= 2A^[P++(g) + 7^P-p-(g)]. (8.5) 

Using the expression of Cppi (Appendix ^ , we obtain 



ImTrqX; 



27r 



yy glN{Q) 



An^N 



InI- 



N-i-N 



-4r 



S{uj - VFqx) 

Vpqx 

2 ^2+^2' 



r 



(l-r2)2 (l-r2)2 |7^| 
5{lo-lox) 3(1 -hr2) I-72I 



LUl 



2(l-r2)2 |7( 



(8.6) 



for w, (/a; > and qy = 0. Both spin-wave modes con- 
tribute to the spectral function. The spectral weight car- 
ried by the Goldstone mode diverges as 1/qx as expected 
for a quantum antiferromagnet.Ej 

Eqs. ( 6.23| ) and ( ^.6|) predict that all the spectral 
weight is carried by the in-phase modes, i.e. the gapped 
sliding mode and the gapless spin-wave mode, whenever 
both SDW's have the same amplitude (7 = C = il)- 
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IX. HELICOIDAL PHASE 

The analysis of the hehcoidal phase turns out to be 
much simpler than the one of the sinusoidal phase. In 
the helicoidal phase, the mean-field gap does not depend 
on the transverse momentum ky, which significantly sim- 
plifies the computation of the collective modes. In this 
section, we shall describe the properties of the helicoidal 
phase, but skipping most of the technical details of the 
derivation. 



where the last line is obtained using the gap equation 
9A_E/9A* -|- = 0. From the latter, we also deduce (writ- 
ing A„, = |A„<,|e**=-) 



V'+T - V-] 



' if InIn > 0, 

TT if iNlff < 0, 



and 



<P+T = f-]-- 



(9.6) 



(9.7) 



A. Mean-field theory 

The helicoidal phase is characterized by the order 
parameters 



A„,(r) = (0„,(r,r))=A„,e"'^: 



p(a,cr)JVl" 



with 



p{+A)=p{-,i) = +, p(+,i)=M-,T) 



(9.1) 



(9.2) 



In the notations of the preceding sections, this corre- 
sponds to A^^^ = A+i ^ A*_^, 



AN) 



AN) 



0, and A 



{N) 



AN) 



A 



0. The fact 



that some order parameters vanish in the helicoidal phase 
makes the computation of the collective modeSptauch sim- 
pler. The spin-density modulation is given bycd 

(5,(r)) =2|A+||cos(-QAr-r-(^+l) 
+2|A_||cos(Q^-r-(^_l), 

{Sy{r)) =2|A+||sin(-Qjvr-(^+l) 
+2|A_||sin(QjvT-^_l), 

(5.(r)) = 0, (9.3) 

which corresponds to two helicoidal SDW's with opposite 
chiralities. The mean-field action is still given by (3.6), 
but the pairing amplitudes are given by (in the QLA) 



AQCT(k, k') = Ski^k^aQo^aa{ky), 

A (k ) — A piap{a,a)N{kyb~TT/2) 

Aqct = Ip{a,ij)N{92^aa + 93^aa)- 

The ground-state energy reads 

-'-nl rv 



(9.4) 



AE ^ 



^p{a,^)N 

N{0)y 



^ + |A.Tpln^ 



2Eo 

AaT 



i7V(0)VlA, 






qTI 



(9.5) 



Without loss of generality, we can choose the order pa- 
rameters real, and introduce 



A_t 



A+T 



^=A 



(9.8) 



+T 



7 and 7 satisfy Eqs. ( 3.2C| ). 

The mean-field propagators are then 



GaaiK^) = 



FaaiK^) 



'ILO - eaa{kx) 



C.2 + 4^(fc,) + |A„.|2' 
A ioip(a,a)N(kyb-T!/2) 



^a(j\kx 



|A„J2 



(9.9) 



211/2 



The excitation energy Eaa{^) = [(^^.A^x) + |A„o-P] 
does not depend on ky . The spectrum is shown in Fig. ||. 
Note that all the branches are gapped, due to the pres- 
ence of two SDW's. In contrast, in presence of a single 
helicoidal SDW, some branches would remain gapless. 

As in the sinusoidal phase, we can compute the mean- 
field susceptibilities. They are given in Appendix y. 



B. Collective modes 



Considering phase fluctuations only, we have 

77„,(r,r) = A«,e'"Q-<°-)""-+^'^°-("-'^) - A„,(r,r) 

~iA„^(r)v3„<,(r,T) (9.10) 



to lowest order in phase fluctuations. Following Sec. IV, 
we then obtain the effective action 



q \ ^ 

(9.11) 



V- 



^('?) = 252A^t(',( 



7r(l -C++ -c__) 



c+H 



2__) 7^(1 -c__ 



(9.12) 
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The dispersion of the collective modes is obtained from 
detV~^{q) = 0. We find a Goldstone mode satisfying 



VFQx, <^+T ^^-V 



(9.13) 



This mode corresponds to a uniform spin rotation. 
Equivalently, it can also be seen as-a translation of the 
two SDW's in opposite directions.c^l Thus it combines 
characteristics of the two Goldstone modes of the sinu- 
soidal phase. 

There is also a gapped mode with dispersion law u'^ — 



9 9 



^2 I 
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^ 52 1- 
The oscillations satisfy 



A+|A_-r 



In In 



(9.14) 






7^(1 - r^y + 2r7(7 + ry + 2r^{l + -fry 



72(1 - r2)2 - 2(7 + r)2 - 2r2(l + -fry 



(9.15) 



When the two SDW's have the same amplitude (7 = 1), 



ip+^/ip-] 



-1. As for the sinusoidal phase, the gapped 



mode is found to lie above the quasi-particle excitation 
gap. 

The similarities with phase modes in two-band or bi- 
layer superconductors are more pronounced in the he- 
licoi dal p hase than in the sinusoidal phase. In fact, 
Eq. ( 9.14 ) gives exactly the gap of the phase modes in 
these superconducting systems. I n th e helicoidal phase, 
the order parameter Aaaiic) [Eq. (9T)] and the effective 
potential A„<,(r) = 92Aaa{r) + gse'^^'^^'^Asair) have 
only one component '^ g'^'^p*"")"''. Thus, the helicoidal 
phase is more BCS-like than the sinusoidal phase. This 
explains why the quasi-particle excitation spectrum does 
not depend on ky, and also the absence of a 7 -depe ndent 
factor in the energy of the gapped mode [Eq. ( |9.14 )]. We 
have shown in Sec. VI that in the sinusoidal phase only 



the largest gap remains tied to the Fermi level in presence 
of condensate fluctuations. This yields a modification of 
the condensation energy, which is at the origin of the 
7-dependent factors in the expressions of uiq and toi . 



C. Spectral functions 

In the helicoidal phase, it is not clear whether the 
collectiiie modes should be seen as sliding or spin-wave 
modes.L3 In this section, we compute the spectral func- 
tions Re[CT(ti>)] and ImTrqx"^"' to answer this question. 



1. Optical conductivity 
The charge and current operators are given by 



Pdw 



(i",t) = — a:r(^^+t(r,T) -(^_|(r,r) 



(9.16) 



JDw(r,r) = -^dr(0+'f{r,T) -^_|(r,r)j, (9.17) 
where 



~ . X g2Aacr'Paa{T^,T)+g3^aaVaair,T) 

ipaa{r,T) = (9.18) 

g2Aaa + gs^aa 

is the fluc tuati ng phase of the effective potential A^a- 
From Eq. (9.17) we can calculate the current-current cor- 
relation function, which yields the conductivity 



Re[a{uj) 



'-S{uj ± UJ2)- 



(9.19) 



The Goldstone mode does not contribute to the conduc- 
tivity and is therefore a pure spin- wave mode (in the limit 
q -^ 0) . This result could have been anticipated since the 
condition i^_|_| = (p _| im plies (^+f = if-i and the vanish- 
ing of the current (9.17). Thus all the spectral weight is 
pushed up above the quasi-particle excitation gap at fre- 
quencies of the order of 072 . This implies that there is no 
low-energy collective mode corresponding to a uniform 
sliding of the condensate. Since non-linear conduction in 
SDW systems results from the depinning of such a mode 
(from the impurity potential) above a threshold electric 
field, we conclude to the absence of non-linear conduc- 
tivity in the helicoidal phase. 

In Ref. ^ we have shown that the helicoidal phase 
is characterized by a vanishing QHE {(Jxy — 0) and a 
kinetic magneto-electric effect. By studying the collec- 
tive modes, we have found a third possibility for the ex- 
perimental detection of the helicoidal phase, namely the 
absence of non-linear conduction. 



2. Spin-spin correlation function 

We consider the spin-spin correlation function 

X^.(r,r;r',T') = (S^r, r)5,(r',T')). (9.20) 



As in Sec. VIII , we shall compute the partial trace 
TrqX^i/ , where the momentum q of the collec tive mode 
is held fixed. [Again, we drop the last term of (3^).] We 
find 



a 
TrqXa;^ = 0. 

This yields the spectral function 



(9.21) 
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ImTrqX; 



rot 


■K 


A+TA_t 


fM 


2glN{Q) 
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5{uj - W2) 



072 



4r 



(1-r 




7' + C' 



(l-r2)2 (l-r2)2 7C 



(9.22) 



for w, ga; > and qx — *■ 0. Both modes contribute to 
the spectral function. Although the Goldstone mode is 
a pure spin-wave mode, the gapped mode has character- 
istics of both a spin wave and a sliding mode, as can be 
seen from the spectral function. 



X. CONCLUSION 

We have studied the long-wavelength collective modes 
in the FISDW phases of quasi- ID conductors, focusing 
on phases that exhibit a sign reversal of the QHE (Rib- 
ault anomaly). We have recently proposed that two 
SDW's, with wave vectors Qat = (2/c_f -|- NG, Qy) and 
Q^ = {2kp — NG, —Qy), coexist in the Ribault phase, as 
a result of umklapp scattering. When the latter is strong 
enough, the two SDW's become circularly polarized (he- 
licoidal SDW's). The presence of two SDW's gives rise 
to a rich structure of collective modes, which strongly 
depends on the polarization (linear or circular) of the 
SDW's. 

Regarding the sliding modes, we find that the out- 
of-phase oscillations are gapless in the long-wavelength 
limit. The fact that this Goldstone mode corresponds 
to out-of-phase (and not in-phase) oscillations is related 
to the pinning by the lattice (due to umklapp processes) 
that would occur for a single commensurate SDW. The 
other sliding mode is gapped and corresponds to in-phase 
oscillations. In Bechgaard salts, this mode is expected 
to lie above the quasi-particle excitation gap and should 
therefore be strongly damped due to the coupling with 
the quasi-particle excitations. In the helicoidal phase, 
there is no low-energy sliding mode, since the Goldstone 
mode is a pure spin-wave mode. [For a helicoidal SDW, 
one cannot distinguish between a uniform spin rotation 
and a global translation, so that we cannot classify the 
modes in sliding and spin- wave modes.] 

The dissipative part of the conductivity, Re[cr(ijj)], ex- 
hibits two peaks: a low-energy peak corresponding to the 
Goldstone mode, and a (broader) peak at high energy due 
to the incoherent gapped mode. The low-energy spectral 
weight is directly related to the ratio of the amplitudes of 
the two SDW's that coexist in the Ribault phase. When 
the umklapp scattering strength ((73) increases (experi- 
mentally this corresponds to a pressure decrease), spec- 
tral weight is transfered from the low-energy peak to high 
energies. Above a critical value of 33, the sinusoidal phase 
becomes unstable with respect to the formation of a he- 



licoidal phase. At the transition, the low-energy spec- 
tral weight suddenly drops to zero (Fig. ^), since there 
is no low-energy optical spectral weight in the helicoidal 
phase. Thus, the formation of the helicoidal phase can 
be detected by measuring the low-energy optical spectral 
weight. We also note that the absence of a low-energy 
sliding mode_.means that there is no infinite Frohlich 
conductivityciJ in this helicoidal phase, which is there- 
fore a true insulating phase even in an ideal system {i.e. 
with no impurities). In a real system (with impurities), 
this implies the absence of non-linear dc conductivity. 

The spin-wave modes exhibit a similar structure. The 
in-phase oscillations of the two SDW's are gapless (Gold- 
stone mode), while the out-of-phase oscillations are 
gapped. In the helicoidal phase, the Goldstone mode is a 
pure spin-wave mode, but the gapped mode contributes 
both to the conductivity and to the transverse spin-spin 
correlation function. 

As discussed in the Introduction, these conclusions rely 
on a simple Fermi surface (Fig. Ilia), which does not neces- 
sarily provide a good approximation to the actual Fermi 
surface of the Bechgaard salts. They should therefore 
be taken with caution regarding their relevance to the 
organic conductors of the Bechgaard salts family. How- 
ever, we have derived a number of experimental conse- 
quences that should allow to test our theory. In the si- 
nusoidal phase, we predict a possible reentrance of the 
phase iV = within the cascade. EJ The low-energy peak 
in the optical conductivity Re[(T(a;)] carries only a frac- 
tion of the total spectral weight uj^/A. This fraction 
should decrease with pressure. At low pressure, the si- 
nusoidal phase may become helicoidal. The, helicoidal 
phase is characterized hy a vanishing QHEJ^J a kinetic 
magneto-electric effect,EJ and the absence of low-energy 
spectral weight in the optical conductivity as well as the 
absence of non-linear dc conductivity. In the alterrLaaj 
five scenario proposed by Zanchi and Montambaux,E3 
the Ribault phase does not exhibit any special features 
compared to the positive phases, aparjti'rom the unusual 
behavior of the magneto-roton modes,c3 but these modes 
have not been observed yet. 

We expect our conclusions regarding the structure of 
the collective modes and the associated spectral func- 
tions to hold for generic two-SDW systems and not only 
for the FISDW phases that exhibit the Ribault anomaly. 
It is clear that the existence of four long-wavelength col- 
lective modes (two spin-wave and two sliding modes) re- 
sults from the presence of a second SDW, which doubles 
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the number of degrees of freedom. Two of these modes 
should be gapless as expected from Goldstone theorem in 
a system where two continuous symmetries are sponta- 
neously broken: the translation symmetry in real space 
and the rotation symmetry in spin space. 

This belief is supported by the striking analogy with 
collective modes in other systems—like two-band, bi- 
layer, and d + id' superconductorsjlJTJ and, to a lesser 
extent, plaanon modes in semiconductor double- well 
structures.Ej While phase modes are in general difhcult 
to observe in superconductors, since they do not directly 
couple to external probes (see however Ref. |l^), collec- 
tive modes in SDW systems directly show up in response 
functions like for instance the dc and optical conductivi- 
ties. 

Finally, we point out that we expect a similar struc- 
ture of collective modes in the FISDW phases of the 



J 



compound (TMTSF)2C104. Due to anion ordering, the 
unit cell contains two sites, which leads to two electronic 
bands at the Fermi level. This doubles the number of 
degrees of freedom and therefore the number of collec- 
tive modes (without considering the possible formation 
of a second SDW due to umklapp scattering). As for 
phonon modes in a crystal with two molecules per unit 
sites, we expect an acoustic (Goldstone) mode and an 
optical mode. 
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APPENDIX A: MEAN-FIELD SUSCEPTIBILITIES x AND x 



Using the expression (3.21) of the mean-field propagators G and F, we obtain 



T 



Xa<7,aa{<i+aQpN,q' +aQp'N,(^,y) = ~'^q,q' r r y^, ^ Gqct (k, tj)G'aa (k - QfQo - q + QfTlG, LJ - LU^) 



x-^y 



li,u n — — C30 



K,LJ n — — OQ 

xIpN+n{n/b + <z,)Vjv-„(vr/6 + qy)e-^^{v+p')N{kyb-,yb/2-./2) ^ (A2) 

Here we have assumed that \qx\^ Wx\ ^ G, which is the case of interest when studying the low-energy fluctuations 
around the mean-field solution. 

In the QLA, we retain only the term n = in the above equBjtipns, since the other terms are strongly suppressed 
when ujc S> T, |ApAr| in the limit of long- wavelength fluctuations Ej Restricting ourselves to qy — Q, we then obtain 



Xacr,acr(q + "QpTV, Q + "Qp'TV, ^i/) 



T 



Ajf J-J 



x-^y 



■ J2 Gaa (k, Lo)G^s (k - aQo - q, c^ - c^.)/pjv/p'Are-*"(P-P')^('=«''-^/2) , 



(A3) 



Xacr,aiy{<l + aQpN , q - aQp' N , ^i^) 



T 



LxJ-ii 



■ J2 Fo.a (k, C^)^a.(k - q, C^ - C^,)/pAr Vjve-'"(P+P' 



)Ar(fc„6-7r/2) 



(A4) 



Performing the sum over kx, we obtain at T = and in the limit \vFqx\^ \^v\ ^ |AQCT(^y)| 



Xa<y,aa{(l+ aQpW, Q + aQp'W , W^) = ^p7V-^j 



p'N 



XQ<T,a5-(q + aQpW. q - aQp' TV, ^v) = IpNip' 



N 



N{0) 

2N\ 



^e- 



(p-p')N(kyb-TT/2) 



In- 



2E. 



2 1 „,2 „2 



1 ujt + Vpq: 



E 



^-ia{p+p')N{kyb-7r/2) 



\A„a{ky)\ 2 6|A„,(fcy)|2 

1 vj,ql+ujl 

2 12|A„,(fc^)I2 



, (A5) 
(A6) 



Eqs. ( |A5| ) and ( [A^ ) are obtained by expanding to flrst order in u;l/\Aaaiky)\'^ and Vpql/\A.aa{ky)\'^ . This calculation 
is standard w hen ev aluating the long-wavelength collective modes of a SDW system and can be found for instance in 
Ref. E3. Eqs. (3.29) are then obtained by summing over ky. 
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It is useful to introduce the notations 

Cpp' =52[x+T,+T(q + QpW:q + Qp'N^L^u) -x+T.-i(q + Qpw.q- Qp'n,l^u)] = cpp> + 6cpp>, (A?) 



where Cppi = Cpp'|q=i^^=o and Jcpp' are deduced from Eqs. ( 3.29| ). 



The gap equation can be written as a function of the static (w^ — 0) susceptibihties x (or c). From AQ,cr(r) 
{Oaaij)) = TJ2ujPacr{r,r,u!), wc deducc 



ApAT = /pAT - 



T 



E 



-i(p-p')Af(fca''-'r/2) 



This can be rewritten in terms of the mean-field propagators 



-A„/ 



p'N- 



(A8) 



^pAf — —IpN 



T 



LrLoi 



E' 



-i(p-p')A'(fe«6-^/2) 



[G+t(k,c.)G_x(k-Qo,c^)-|F+T(k,c.)|2]Ap, 



AT 



7 , 7 \x+1, + T{QpN,Clp'N,^iy — 0) - X+h-i{QpN,-Qp'N,^i> — 0) 

^-^ i-n'AT L 



T -'p'JV 



Using the relation between Apjv and Ap^r, we obtain 



A„, 



p'N- 



(A9) 



^7V 



rA 



TV 



c++An 



A 



AT 



rA 



TV 



c__A 



AT 



(1 - r2)/A. (1 - r2)/^ 



Af 



^JV 



(1 - r2)/^ (1 - r^)l 



N 



In 



From (AlO), we deduce the relations 



C++ 



1 + 7r 



c„_ 



7 + r 



which gives, by eliminating 7, 



(1 - c++)(l - c — ) - r^c++c- 



0. 



(AlO) 



(All) 



(A12) 



Eq. (A12) is nothing but the gap equation rewritten in terms of the static mean- field susceptibilities 



In the study of the collective modes, the natural quantity to consider is not x but the susceptibility x defined by 



Xaa,aV'(r,T;r',r') = (0„,(r, T)0;,,,(r', r')) - (0„,(r,r))(d:,,,(r',T')), 
where Oa{r,T) = g{r)Oa{r,T). x is related to x by 

Xa<T,Q'a'(q, q',Wj.) =g|xacr,a'o-'(q,q',w^) + gixa<T,a' a' (q ~ aAkp , q' - a'4kF,UJ^) 

+g2g3[Xac7.a'a' (q - a4fc_F, q', UJ„) + Xaa.a'a' (q, q' " a'-ikp, l^y)] 

and satisfies 

XaT,a'T(q + "QpJV,q + a'Qp'JV,Wiy) - XQT,a'-l(q + "Q!3JV,q- a'Qp'JV,^!.) = g2Sa,a'{Cpp' +r'^Cpp') 

iSS^a^a' \Cpp' ' ^pp' ) • 



(A13) 



(AM) 



(A15) 



Eq. ( [A15| ) is used in Sec. |lv 



APPENDIX B: CHARGE OPERATOR pDw 



Using 



C^(r,T;r',r') 



G„T(r,r;r',T') ^„|(r, r; r',T') 
^a^(r,r;r',T') Ga|(r, r; r', r') / ' 



where G and F are the mean- field propagators (see Sec. IIIB), we rewrite Eq. ( 6.12| ) as 



(Bl) 
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PDw(r,r) = -^ /dV'dT'[G„T(r,T;r^r077„T(r^T0Fax(r^T';r,T)+F„t(r,T;r^r')CT(r^^')GaT(I■^T';r,T) 
+i^ax(r,r;r',T')ryaT(r',T')Gai(r',r';r,T)+Gax(r,r;r',T')CT(i-',T')i"aT(r',r';r,r) . 



If we consider phase fluctuations only, rjaa is given by (iA), which gives using ( |6.8| ) 

p 
Defining the fluctuating phases (faa of the effective potential Aqct by 



A 



A 






pN 



we have 






The relation between (paa__ and (/Jqo- is given by (6.15) to lowest order in phase fluctuations. 
From Eqs. (B2) and (B5), we deduce 



r)+i(q'-r'— oj^r') 



LxLy '—' IpN ■^ J J 

" p '^ a,q' 

X |e-Qp«-r' [G„^ (r, r; r', T')F^i{v' , r'; r, r) + F^^{v, r; r', r')Gai(r', r'; r, r)] 

where q = {q^, qy = 0) and q = {qx.,uj„). Using 

/" dV dr / dV dr' e-''(q-'--'^.^^)+»(q'-r'-^:.r')+,:aQ^„.r'(^^^ ^j.^ ^. r', r')Fax(r', r'; r, t) = 

k.cj 

/" dV dr I d?r' dr' e~^i^■''--^'^^)+^i^'■r'-^'y)+^aQ,^■r' p_^^j.^ ^. ^.r r')Gsiir', r'; r, t) = 

h^'IpN J2 ^ai(k,c.)GaUk - q, w - a;,)e*"P^('^-='^-^/2), 

k,a' 

k.cj 

(5g,g'/pA'^Gaj,(k,w)Fa|(k + aQo -q,t 



,uj - uj^ e 



-iapN{kyb-TT/2) 



k.u 



we obtain 






x< e 



a,p ^ k,cj 

.apAr(fe„6- V2) [G^^ (k, ^)i;^_^ (k - aQo - q, ,, - c^,) + F^x (k, c^)Gax (k - q, c^ - c^.)] 



-iapN{kyb-TT/2) 



(B2) 
(B3) 

(B4) 

(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 



[Fa|(k,w)Ga|(k-q,w - w^) + Gax(k, w)i^c<T(k + aQo -q,cj -tJt.)]|. (Bll) 
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To lowest order in vpq^ and w^, we have (for qy — 0) 



n— > Gaai^, i^)Faa^ - aQo - (\, UJ - LO i,) = -N{Q)—^ -——-{iujy + avpqx), 

hi,.. ■^— ^ p I '^ ' ' ^ I ' 



bL 



kx,UJ 



8|A„T(fc,)|2 



—— y i^ag-(k,cj)G'affk- q,w- Wi,) = -iV(0) ' i-iui, + avFqx)- 



6L 



8|A„t(fcy)|2 



(B12) 
(B13) 



This leads to 



PDw(?) = ^ E ApArat.^ g,e-P^( 



4iV_L 



p^a,ky 






Performing the sum over fc„, we finally obtain 



.-s . qx / ..{pN),~. ..(pN),-. 



which yields Eq. ( 6.13| ) 



(B14) 



(B15) 



APPENDIX C: MEAN-FIELD SUSCEPTIBILITIES IN THE HELICOIDAL PHASE 



In the helicoidal phase, the mean-field susceptibilities are given by 

Xaa,a(7{(l + C^Qp(a,o-)Af 7 Q + OiQp{a,a)N t ^i^) — Ip{a,a)N — ^ — 

-'-p{a,ry)N n 



ln4^° 



Xaa,aa{(l + aQp(Q,CT)7Vj Q " aQp(a,cr)Af j ^y 



1 ml + vlql 
,Aa.| 2 6AL 

1 cj^_+v|g| 

2 12A2 



(CI) 
(C2) 



for qy = and jwi/j, Ui^jq^:! -C |Aq,o-|. Introducing the notations 



C++ = g2[x+h+]i^ + QN,(l + QN,uj„) 






the gap equation reads 



(l-c++)(l-c__)-r2c++c__=0, 



where C++ = 1/(1 + 7r) and c — = 7/(7 + r) [cpp = Cpp|^=o)- 



(C3) 



(C4) 
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(a) 



> 
Q2 



(b) 




FIG. 1. (a) Fermi surface deduced from the dispersion law (1.2) with k = 0. 
Qi and Q2 are coupled by umklapp scattering, (b) Fermi surface for k — n/4. 
By umklapp scattering, it couples to Q2, which is not a good nesting vector. 



At half-filling, the two 'best' nesting vectors 
There is only one 'best' nesting vector (Qi). 
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FIG. 2. Transition temperature T^ between the metallic phase and the FISDW phases in presence of umklapp processes. 
The vertical lines are only guides for the eyes and do not necessarily correspond to the actual first-order transition lines between 
FISDW phases, (a) P3/32 = 0.03. [(73 and p2 are the strengths of normal and umklapp processes, respectively (see Sec. H).] 
The shaded area corresponds to the Ribault phase N = —2. (b) gs,/ g2 = 0.06. Two negative phases, N = —2 and N = —4, 
are observed (shaded areas). The phase TV = —2 splits into two subphases: helicoidal (dark shaded area) and sinusoidal (light 
shaded area). 




-0.5 



0.5 



FIG. 3. Quasi-particle excitation gap 2jA-|-|(fcj;)| in the phase A'^ — —2 for I7I — 1/3. 



23 



tJ/2|AJ^ 
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"^ = ^>Wx + "O 



«Ffe 



FIG. 4. Dispersion laws of the two sliding modes in the Ribault phase. The gapped mode is generally expected to lie above 
the quasi-particle excitation gap. The latter varies between 2|A]v|(l — I7I) and 2|A]v|(l + I7I) (Fig. H). 
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FIG. 5. Fraction of the total optical spectral weight ijJp/4, carried by the low-energy collective modes. In the sinusoidal phase 



(I7I < 7c ~ 0.4), the low-energy sliding mode carries the fraction 3(1 
helicoidal phase, there is no spectral weight at low energy (see Sec. IX). 



7 )/(3 + 57 ) of the total spectral weight. In the 



0.5 



0.03 0.05 0.07 0.09 

J" = 53/52 

FIG. 6. I7I (dashed line) and low-energy optical spectral weight (solid line) as a function of r = 93/92 
{H = 10 T). The transition from the sinusoidal to the helicoidal phase occurs for r ~ 0.06. 



in the phase N — —2 
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*" f^x 



FIG. 7. Quasi-particle excitation spectrum in the helicoidal phase. The sohd (dashed) hne corresponds to up (down) spins. 
For clarity, we have not shown the Zeeman splitting. 
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